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ABSTRACT 


This  report  deals  with  the  theory  and  design  of  antennas  printed 
on  or  embedded  in  a  grounded  substrate.  A  theoretical  approach  is 
implemented  which  accounts  accurately  for  the  physical  effects  involved 
including  surface  waves.  The  Green's  function  has  been  obtained  by 
synthesizing  the  fields  of  Hertzian  dipoles  which  are  oriented  in 
arbitrary  directions  and  which  are  printed  on  or  embedded  in  the  sub¬ 
strate.  By  considering  Pock! ington ' s  integral  equation  with  proper 
choice  of  expansion  and  testing  functions,  solution  for  current  dis¬ 
tribution  and  input  impedance  is  obtained  by  matrix  inversion. 

As  an  example,  thin-type  printed  and  embedded  circuit  antennas 
with  symmetric  or  asymmetric  excitation  are  considered.  Antenna 
currents,  input  impedance,  bandwidth  and  resonant  resistance  are 
obtained  for  a  variety  of  antenna  arrangements.  A  serious  amount  of 
effort  is  also  being  placed  in  evaluating  the  importance  of  higher 
order  surface  wave  modes  which  are  determined  by  the  relative  dielec¬ 
tric  constant  and  the  thickness  of  the  substrate. 
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CHAPTER  1 


introduction: 

Monolithic  or  hybrid  integrated  circuits  are  finding  increased  use 
in  the  microwave,  millimeter  and  far  infrared  frequency  ranges.  There¬ 
fore  the  development  of  antennas  which  are  amenable  to  integration  with 
other  printed  circuit  elements  is  of  significant  technolog ical  impor¬ 
tance  . 

Printed  circuit  (microstrip)  antennas  were  apparently  introduced 
first  in  the  early  1950's  [1],  [2],  However,  their  inherent  advantages 
(conformality  to  a  given  surface,  light  weight,  negligible  volume,  in¬ 
expensiveness),  were  not  put  to  widespread  practice  until  the  1970' s 
[3]  -  [20],  The  need  to  integrate  microwave  components,  including 
antennas,  necessitated  the  development  of  substrate  materials  with  the 
appropriate  thermal,  mechanical  and  electrical  properties  in  the 
desired  frequency  ranges  [21]  -  [23]  .  The  environmental  and  technolo¬ 
gical  constraints  having  been  resolved,  the  task  remained  to  develop 
analytical  methods  which  would  provide  accurate  desicn  criteria. 
Extensive  efforts  have  been  expended  to  provide  models  which  predict 
correctly  the  electrical  characteristics  of  microstrip  antennas.  These 
methods  rely  either  on  a  transmission  line  model  of  the  printed  circuit 
antenna  [3]  -  [7]  or  an  open  resonator  model  [f]  -  [22'].  The  former 
approach  gives  a  heuristic  explanation  of  the  radiation  properties  of 
the  antenna  while  the  latter  provides  a  more  accurate  prediction  of  tie 
antenna  charact  er  ist  ic  s  .  Both  models  applp  main.lp  to  the  dominant 
resonator  mode  and  their  accuracy  is  quest  ion, a-' 1  e  for  highet  order 
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•  ,  tspiiially  because  tiny  do  n,  !  account  for  the-  excitation  of 
surface  waves. 

Surface  waves  are  important  contributors  to  the  printed  circuit 
antenna  current  distribution  as  well  as  input  impedance  characteristics 
(24],  In  addition,  since  surface  inodes  decay  as  the  inverse  square  root 
of  distance  from  their  source  they  can  be  significant  to  mutual  imped¬ 
ance  computation  [25],  a  calculation  of  importance  in  phased  array  de¬ 
sign.  It  has  been  established  that  regardless  of  how  thin  the  substrate 
is,  the  dominant  surface  wave  mode  is  always  excited.  The  efficiency  of 
launching  of  this  mode  depends  on  the  thickness  of  the  substrate  and  its 
dielectric  constant.  As  more  energy  is  trapped  in  the  substrate  the 
microstrip  antenna  becomes  less  efficient.  In  many  applications,  such 
as  ir.  the  millimeter  or  far  infrared  region  [2?]-  [3C],  today's  technol¬ 
ogy  provides  substrates  which  are  several  wavelengths  thick.  This  per¬ 
mits  many  surface  modes  (?E  and  TM">  tc  exist  in  the  substrate,  further 
complicating  the  design.  Tnese  modes,  too,  can  cause  impairment  of 
ef  f  ic  ien.cy . 

It  becomes  evident  fron  this  discussion  that  a  theoretical 
approach  must  he  implemented  which  accounts  accurately  frr  th«.  physical 
effects  involved,  including  surface  waves.  Such  ar  aprronr:  exclude.- 
either  of  the  previously  mentioned  techniques  and  relies  on  treating 
the  microstrip  element  as  an.  antenna  rather  than  as  a  transmission 
line  section  or  as  a  resonator.  This  requires  that  the  antenna  current 
distributin'  he  obtain*  first  by  solver,  a  t  wo-d  ir.ensionr.  1  •  ock  ling  tor. 
integral  equation.  The  C-reer's  function  m  til'  case  can  be  obtained 
b;-  s  ynt  f.es  ir  in.  the  fields -of  Fertrinr  d.rl-<  whi:  1  art  oriented  ir 


arbitrary  direction'  and  whirl  arc  printed  or.  the  substrate,  thus 


account  i nr,  properly  for  all  the  boundary  conditions  pertinent  to  the 
problem . 

An  analytical  solution  of  the  two  dimensional  Pockl ington ' s  inte¬ 
gral  equation  is  precluded  due  to  the  immense  complexity  of  the  problem. 
One  must  resort  to  numerical  techniques,  with  the  integral  equation 
discretized  and  the  current  distribution  obtained  by  matrix  inversion. 

A  numerical  method  which  has  found  widespread  and  successful  use  for 
the  solution  of  Pocklington  type  integral  equations  is  the  method  of 
moments  (24]  -  [27].  It  has  the  disadvantage,  like  other  numerical 
techniques,  of  requiring  extensive  computing  time  in  some  applications, 
e.g.,  in  the  computation  of  electromagnetic  scattering  from  objects 
large  compared  to  wavelength.  For  the  present  application,  the  Green's 
function  pertinent  to  the  problem  is  given  by  Sommerf eld-type  integrals 
which  require  special  integration  techniques  when  field  and  source 
points  are  both  on  the  substrate  [24  ] ,  [25],  [31].  At  present, 

evaluation  of  these  integrals  is  numerically  time-consuming,  thus 
making  the  applicability  of  the  method  of  moments  to  arbitrarily 
shaped  planar  printed  circuit  antennas  difficult. 

In  the  present  work,  thin  wire-type  printed  circuit  antennas  with 
symmetric  or  asymmetric  excitation  are  considered.  The  excitation  is 
assumed  to  be  caused  by  a  unit  voltage  generator  and  the  printed  vires 
are  taken  to  be  thin  compared  to  wavelength.  The  azimuthal  asymmetry 
of  the  current  distribution  is  assumed  tc  be  a  second  order  effect. 

The  antenna  currents,  input  impedance,  bandwidth  and  resonant  resistance 
are  obtained  for  a  variety  of  antenna  arrangements.  A  serious  amount 
of  effort  is  alsr  being  placed  in  determining  t  he-  importance  of  higher- 


3 


ori’ti-  surface  wave  modes  so  that  a  ecu  net  trade-off  analysis  of  band¬ 
width  with  substrate  thickness  can  be  implemented. 


CHAPTER  2 


GREEN'S  FUNCTION  FOR  PRINTED/EMBEDDED 
DIPOLES  AND  POCKLINGTON'S  INTEGRAL  EQUATION 


2-1.  DERIVATION  OF  GREEN'S  FUNCTION 

This  chapter  presents  the  development  of  the  Green's  Function  per¬ 
tinent  to  the  problem  of  linear  antennas  printed  on  or  embedded  in  a 
grounded  dielectric  substrate  of  thickness  b  and  relative  dielectric 

constant  e  . 

r 

In  order  to  formulate  the  Green's  function,  an  elementary  hori¬ 
zontal  electric  dipole  (HED)  is  considered  to  be  at  (x',  y',  -h)  as 
shown  in  Figure  2.1.  The  assumed  time  dependence  is  e^ u;t  and  it  is 
suppressed  throughout  the  report.  Maxwell's  equations  take  now  the 


following  form  in  region  II  (z  <  o) 

ri  rtd  -id  .  ■*  d 

»  x  h  *  J  +  jwr  e  I 


(2.1) 

(2.2) 


*  E  =0 


(2.3) 


(2.4) 


where  the  superscript  d  indicates  field  and  source  quantities  ir.  th<. 
dielectric  substrate.  The  following  potential  function  is  now  intro¬ 


duced  , 


rtd  kS  -  d 

-•  h  =  — —  .  x  C 
«  J- 


(2.5) 


where  C-  will  turn  out  to  be  the  Green's  function  ir.  region  II.  Since 
V  ’  x  C-  ^  )  *  0,  equation  (2.5)  indicates  that  (2.3)  is  satisfied.  A 
substitution  of  (2.5)  into  the  cur)  equation  for  E ^  gives  the  following 


ri-sul 1 


V  x  (Ed-k2Gd)  -  0 


(2.6) 


where  k2  =  w2p  e  c  .  Since  V  x  (V4d)  .  0  it  follows  that 
o  r  o 


:d-k2  <*d-  v;d 


(2.7) 


with  4  the  scalar  potential  function  for  region  II.  The  relationship 
between  4^  and  G d can  be  obtained  by  a  substitution  of  (2.5)  and  (2.7) 


into  (2.1),  i.e. 


-*  -*  -*  d  -*■  -*•  d  .  ± 

V  x  V  x  C  -  V  (9  •  G  )  -  V 


2  d 


G  = 


°  Jd  + 


+  k2  G  d  -  Vi  d 


(2.6) 


4  (J  4 

Since  the  curl  of  G  has  been  defined  and,  in  addition,  since  G  is 


arbitrary,  one  can  choose  its  divergence  to  be 
V  •  Gd=  -  <d 


(2.9) 


which  is  the  Lorentz  condition.  The  wave  equation  for  G  can  now  be 


obtained  from  (2.8)  and  (2.9)  in  the  form 
(V‘  +  k-)  Gd  -  — =£ 


(2.10) 


-4  C 

The  Lorentz  condition  enables  one  to  write  in  region  II  for  E  and  K 
the  result 


ant 


E  d  =  V('  •  Gd)  +  k2  C  d 
0 

ud  ’*'■"* 

h  = . v  X  G 

J^o 


c.i:' 

C  .12) 


In  a  similar  manner  the  electromagnetic  field  in  region  I  is  given  by 


E  ■=  V("  .  t)  +  1  2 


(2.131 


an: 


K  *  -  n-1—  T  x  C 
>4 


(2.14) 


where  is  the  free  space  wavenumber,  i.e.  /  ~ 


c  .  G  is  thi 
c  o 
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C?  1 1  r. 1  v  fum  t  ion  for  region  1  am1  it  satisfies  the  wave  equation 

(V~  +  kJ)  G  =  — J  (2.15) 

0  k2 

o 

4 

for  a  given  current  source  J. 

The  solution  which  satisfies  the  wave  equation  and  the  appropriat 
boundary  conditions,  in  each  region  of  interest,  consists  of  two  parts 
the  secondary  solution  which  is  the  solution  to  the  corresponding  homo 
geneous  equation,  and  the  primary  solution  which  is  the  particular 
solution  to  the  wave  equation.  The  complete  solution  is  therefore 
C  =  Gp  +  Gs  (region  I)  (2.16) 


anc 


G  d  (region  II) 


(2.17) 


The  existence  of  the  primary  solution  is  strictly  related  to  the  pre¬ 
sence  of  current  sources  in  the  region  of  interest. 

For  the  cases  considered  here  the  primary  solution  has  an  x-co r- 
ponent  or.lv  [32]-[3-]  giver,  by  (see  Appendix  A'' 


.  d  jL"c  f  y  \  ~u  ’  2  +1 . '  )  d 


where 


and 


o  1  o  n  1 

VJ  \  u  =  [  >  -  V.  I 


:  =  [(>:  -  x'l2  +  (y  -  y  ’ ) 2  ]  '* 


cos:  = 


alser.rc  cf  source.-  ir.  region  7  th»  secondary  components  an 


Due  tc  the 


(2.19) 


5  -  -  -  °  I  F(X)  j  (>p)  e  -Vd> 

SX  4rk2  Jo  0 

o 


C  * - r-  COS  C 

SZ  4n  k2 
o 


f}:l)  J1 


(>c )  e“Vd> 


(2.20) 


f  |  Jo()p)[Fd(X)eUZ4  Rd(X)e'uz]  dX 


M  /*  OP 

f  /  VXp” 

k  ./  o 

Vo  /“ 

—  COS*  I  J1 

k- 


(Xc)  [«d(X)euz+  « d(X)e_UZ]d; 


(2.21) 


(2.22) 


In  the  equations  above  F(X) ,  4(X),  Fd(X),  Rd(X),  4d(X)  and  d(X) 


determined  by  the  boundary  conditions  at  z  -  0,  -b  given  by: 

E  =  Ed,  E  =  Ed.  H  =  Hd,  H  =  Kd  (z  =  0) 
v  v  >:  x  v  y  x  x 


(2.23) 


Ed  =  0  Ed  *  0  (z  =  -b) 
V  x 


(2.24) 


Application  of  these  boundary  conditions  yields 

d  (1  -  cf) 1 2  eUb  sinh(u(b-hX ' 

fjO.b)  •  f 2  ( X  ,b) 

d .  d  .  -2ub 
-  O)  =  4  (/)  e 


(2.25) 


Fd  ( > ) 


Rd(>) 


y  (u-u^)  sinh[u(b-h)] 
u  fjT'.b) 

r  }  _d , . . .  -2ub 

-[  -  e  +  F  O)  ]  e 


(2.2?) 


(2.2F) 


2r ^ >  Sinh [ u (b-h) ] 

f](' 


2'  (1  “t  j.)  ‘  'sinh [u(b-h)  IcosMu)  ) 


flf.b)  5  2(  -  ,V) 


(2 ,3^ 


it'.  Which 


f^(X,b)  =  u^sinhfub)  +  ucosh(ub) 


f^CXjb)  *  Er  u^cosMub)  +  usinb(ub) 


(2.31) 

(2.32) 


With  the  above  expressions  substituted  in  equations  (2.18)  -  (2.22)  the 
final  form  of  the  Green's  function  is 
jwu 


Gv= - °  f  J^(X  p)e_UoZ  slnh[u(b-h) ] 

i  J  o 


x  ™\j0  v  fi(x-b) 


Xd>. 


G  = 
z 


jwv  r 

- ^  ( 1  — c  )cosf  I  J.(>p)e  V  sinh[u(b-h)cosh( 

2rk;  r  J0  fjTxVb)  •  f2"a. 


h)cosh(ub)  .2 

"'TTET 


(2.33) 


(2.34) 


cl-  .-4?  I  JO.)  - 


*  2rk“  '  ° 


L 


sinh[u(z+b) ] • [ucosh(uh)+u  sinh(uh) ] 


S1  6  ,bj 


-  dX,  -b<z<-h 
u  - 


(2.33) 


Gd=  - 
x 


f  ^  ^  sinh  [u(b-h)  ]  [ucosh(uz) -u0sinh(uz)  ]  )  ,  -h'z<C 

2rk 2  Jc  0  ~f1o,T;>  « 


Gd= 

z 


^""o  ,,  x  N  sinh[u(b-h) )cosh[ 

-  ^5  r  /  1  _"“f1“O.Tw;— 


h  [u  ( z+b)  ] .  2  j  . 
b)  '  d' 


(2.36) 

(2.37) 


For  the  case  of  the  HEh  on  the  dielectric  interface  the  Green’s  function 

is  giver,  by  equations  ( 2 . 33  X  ,(2 . 34  >  for  region  (I),  and  by  (1.33''. 

(2.37)  for  region  (II)  with  h  =  0. 

2 -2 .  POCK LI NOTOh ’ ?  I NTEGRAL  EQUATION . 

It  is  assumed  that  the  dipole  shown  in  Figure  2.2  is  very  thin  so 

that  its  radius  y  is  much  smaller  than  tht  vnvelencth  >  in  the  di- 

F 

electric,  i.t. 


,  '■  =  /  /  >  . 
g  P  o  .  r 


(2.36'' 


It 


I:  it  is  assumed  that  the  current  density  on  the-  dipole  is  given  by 


rL/  32Gd  aV  \ 

■  /  /kV  +  . i+  — 

Jo  \  X  9x3*  ' 


Jd  dr ' 
x  x 


It  is  easily  determined  now  that  the  x-component  of  the  electric 
for  the  case  h  *  0  is  as  follows: 


E*(r) 


,(>c  ♦  %  ♦  ‘-M  , 

Jo  \  OX  3x2  3x3z/  x  x 


with  G  ,  Gd ,  G  ,  Gd  defined  previously. 

X  X  2  z 

A  consideration  of  the  following  relationships 


and 


?G 


02 


JG 


ex 


?G 

c>: 


ci- 

r  V 


leads  to 


idy(7)  = 


■[  l‘-“-  •  5  • 


cd) 


G) 


Jd  dr' 
x  x 


J  cr ' 

X  X 


with  G  ,  G  giver,  by  (see  Appendix  El 


G  = 


2  r  V. 


(1  -c 


-i. 


,,  v  sinh  tu (b-hN , 
o  fjOVbl  u- 


sinh  (_z_  b)  ] 


>  d 


j'-0  f  .  -u  z  sir'.''-..'’''  ccsMu'rl... 

- 2  (rr-1}  /  Jo°p)  e  °  UcV.'  ,r>  •  O-.n  C 

2*Vc  J 


(2  .AS) 

f  ield 

(2.46) 

(2.47) 

(2 . 4?  > 

(2 . 4 ?■) 

C.5C' 

(2.511 

(2 .521 
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CHAPTEI 


FORMVI.ATION  OF  THF.  PRORhF.M  FOP 
ETNlfh  1  CVL_ANALYSI? 


3-1 .  METHOD  OF  MOMENTS 

The  purpose  of  this  chapter  is  to  present  the  basic  principles  of 
the  Method  of  Moments  [35],  [371  as  well  as  its  application,  to  the  spe¬ 
cific  problem  of  the  dipole  embedded  in  a  grounded  dielectric  substrate. 

Throughout  this  chapter,  the  radius  of  the  dipole  is  assumed  tc  be 
much  smaller  than  the  wavelength  in  the  dielectric,  so  that  the  current 
may  be  assumed  to  be  distributed  uniformly  around  the  cylinder.  If  ont- 
observes  the  current  distribution  from  an  observation  point  on  the 
cylinder  axis  then  the  distance  between  the  source  and  the  observat ic-r. 
point  reduces  t 

.  =  [(>:  -  ‘  j  ‘ 

It  is  tc  be  emphasi7ec  that  the  dipcle  or  wire  cylinder  is  net 
assumed  tc  be  infinitesimally  thin  and  therefore  the  singularity  which 
arises  ir  the  computet  ion.  of  reactance  is  overcome  vhe-  the  vir<. 
radiu-  goes  tr  rer: . 

The  unifying  concept  of  the  numerical  treatment  of  input  impe¬ 
dance  and  radiation  problems  is  the  method  of  mo-e-ts.  This  general 
approach  tc  this  kind  cf  rrcbie~«  ic  essert  ir  1 ■  a  reduction  of  the 
associate-'  in*.  err..l  ecus:  tc  a  sys  ter  c  f  !i*e.--.c  algebraic  equation, 

ir  h  unknown-  .  Tot  h  unknown  s  art  usual  1*.  coefficients  in  sore  appro¬ 
priate  expansio:  of  the  current  d is tr ihut  i .  •  . 


i  c, 


i 


The  gi-m-ral  method  of  solution  will  be  discussed  In  the  notation 


of  linear  spaces  and  operators  and  hence  the  specific  problem  studied 
here  will  be  put  into  this  notation.  The  integral  equations  derived  in 
Chapter  2  (2.52  and  2.53)  can  be  put  into  the  form  of  an  operator 
equation  using  the  concept  of  linear  vector  spaces  as  follows: 


L  (J)  -  (E  x) 
op  x 


with  L  given  by 


op 

in  air  and 
L 

op 


ft 

/>•< 


lc  C  +  (C 
0  X  3x2  X 


+  ^(cd 

3x2  x 


C)J  (x-d^) 
Gd)"|  (x-dr^) 


(3.1) 


(3.2) 


(3.3) 


in^he  dielectric  substrate. 

In  equation  (3.1)  (E^x)  Is  8  hnowm:  excitation  function  or  source  and 
is  the  unknown  current  distribution.  In  addition  to  the  above  it  is 
necessary  to  identify  the  inner  product  <  J,  (E^x)  >>  which  is  defined 
to  satisfy  in  Hilbert  space  the  following  relations: 


< 

3,  (E  x) 

V 

n 

A 

M 

X  i 

<-.1 

V 

( 3  .  *  i 

A 

P 

(-.1 

X 

+  e (e  x) 

X 

,  h  >  *  a  <3,  h  >  +  £  <  (Ex),  h  > 

(3.5) 

< 

X 

J*,  J  > 

X 

>0  if  J  *  0 

(3.6) 

-0  if  J  *=  0 


where  a,  £  are  scalars;  j  x  h  =  0  and  *  denotes  a  complex  conjugate. 
The  inner  product  may  be  in  the  fern 

<  J,  (E  x'  >  m  f  J  •  (E  >:1  dr  (?.“) 

x.  J  x 

o 

and  since 


G  (r/r’)  -  G  fr'/r) 


(3 .  £ ) 
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Ul 


(3.9) 


a  2  ,  ,2 

OX  <*X 

then  it  can  be  shewn  easily  from  equations  (3.2)and  (3.3)  that 
*  lop  <3>  •  <Ex*>  >  ■  <  3'  Lop  <E**>  » 


(3.10) 


which  means  that  the  integral  operator  L  is  self-adjoint.  Based  on 

op 

this  fact  and  on  the  physics  of  the  problem  it  is  concluded  (see 
Appendix  C)  that  there  exists  a  unique  solution  to  functional  equation 
(3.1)  and  therefore  the  existence  of  the  inverse  operator  results  as  a 
consequence,  i.e. 

-V  _1  -* 

(3.11) 


3  «  1  _1  (Ex) 

op  X 


The  procedure  for  obtaining  a  solution  in  the  form  of  equation 
(3.11)  can  be  divided  into  four  steps: 

1.  Expand  the  unknown  vector  in  a  series  of  basis  functions.  J  , 

r. 

spanninr  3  ir.  the  domain  of  1 

op 

2.  Determine  a  suitable  inner  product  and  define  a  set  of 
weighting  functions. 

3.  Take  the  inner  products  and  thereby  form  the  matrix  equation. 
Solve  the  matrix  equation  for  the  unknowns  3. 


3 - 2  .  GALEREI^’ S  _ME TKT 

For  electromagnetic  input  impedance  and  radiation  problems  a  spe¬ 
cialization  of  the  general  method  of  moments  is  particularly  convenient 
In  the  first  cr  the  four  steps  above  the  current  distribution  function 
is  expanded  ir.  a  finite  series  of  basis  functions  3^  ( j  «  1 ,  .  .  . ,  h'' 
define:’  ir  the  dor.;. ft  of  L 


< 


(3.12) 


This  leads  tc  the  matrix  equation 


Equation  (3.16)  may  be  solved  by  standard  matrix  inversion  techniques 
for  the  unknown  current  vector  [Jj,  i.e., 

[J]  -  [Z]"1  [V]  [3.17) 

where  [Z]  =  generalized  impedance  matrix 

[V]  =  generalized  voltage  vector 

3-3.  CHOICE  OF  BASIS  FUNCTIONS 


For  the  numerical  solution  of  a  given  input  impedance  problem,  one 

of  the  fundamental  decisions  to  be  made  is  the  choice  of  basis  functions 
— ► 

Jj .  Theoret ically  these  are  infinitely  many  sets  of  basis  functions 
but  practically  there  are  only  a  few  which  provide  an  optimum  solution 
to  the  trade-off  between  accuracy  and  computer  time  for  the  specific 
problem.  However,  it  is  true  that  the  closer  the  basis  functions  J 

n 

resemble  the  actual  current  distribution  on  the  radiator,  the  better 
the  convergence  and  stability  of  the  generalized  impedance  matrix. 

There  are  two  classes  of  basis  functions: 

1.  The_  entire-domain  basis  funct ions .  These  are  defined  and  are 

non-zero  over  the  entire  domain  of  L  except  possiblv  for  a 

op  h  K 

countable  union  of  sets  with  measure  zero.  They  also  satisfy 
the  boundary  conditions  for  the  specific  proble-. 

2.  The  sub-<kxr,ain  basis  functions  which  are  non-zero  functions 
over  parts  of  the  domain.  In  the  case  of  wire  radiators  this 
consists  of  dividing  the  antenna  into  overlapping  or  non-over- 
lapping  sections. 

The  implementation  cf  one  or  the  other  kind  of  basis  functions 
depends  on  the  specific  proller.  In  the  cast  of  arbitrarily  shaped 


radiators  the  problem  can  be  simplified  if  sub-domair.  basis  functions 


an.  used  to  approximate  parts  of  t lie*  radiator  by  N  wire  segments.  Also 
in  the  case  of  certain  geometries  (linear  dipole,  open  or  closed  cir¬ 
cular  loop)  if  the  segments  representing  the  radiator  are  identical 
then  there  is  an  appreciable  reduction  in  computer  time.  This  is  due 
to  the  fact  that  the  formed  integrals  extend  on  the  sub-domains  and  not 
on  the  entire  domain. 

For  the  problem  of  the  printed/embedded  dipole  the  most  appropriate 
set  of  basis  functions  is  found  to  be  of  the  following  form  (see  Figure 


3.1) 

-r-  -v 

J .  =  x 
1 


sin [ k(x ' -x j _j)  ] 
sink.t 

x 


+  P. 
1 


sin  [ k(x -j 4}  -x’) 
sink-' 

x 


6 (z+h) 


(3.16) 


wi  t  h 


P. 

3 


-1 


1  y- .  i  <  x  ’  <  x . 
J-1  j 


and  P .  = 

i 

These  basis  func 
will  be  used  for 
sec  t  ion . 


0  elsewhere 


1  x  .  <  x  ’  <  x  .  , . 
J  1+1 


0  elsewhere 

tior.s,  together  with  the  inner  produc: 
the  formation  of  the  matrix  ecuatior. 


(3.19) 


(3. 2C) 


def ined  earl icr  . 
in  the  next 


3-4.  FORMVLATIPN*  OF  THE  MATRIX  EQUATION 

As  mentioned  previously,  the  choice  of  basis  functions  is  deter¬ 
mined  by  many  factors  dictated  by  the  prohler.  under  consideration.  For 
the  case  of  many  parallel  dipoles  printed  on  or  embedded  in  the  di¬ 
electric,  (see  Figure  3.2)  the  basis  functions  can  be  modified  slightly 


tc  be  of  the  form 


(3.21) 


r 


J  =  f(x,x';j>  6(v'-y  )  d(z'-z  ) 

j  '  p  q 

In  this  case,  the  current  distribution  Inside  a  surface  surrounding  all 
the  dipoles  will  be 

M(p,q)-1 


j  -  E  E  if”  1  pq  Epq 

p  q  j=N(p.q)+i 


(3.22) 


By  substituting  (3.22)into  (3.1),  and  by  using  (3.2)  and  (3.3)  it  car. 
be  found  that  the  electric  field  in  air  and  in  the  dielectric  is  as 
f ol lows : 

•  <  /  \  < 

1 


pq 


M(p, ql -1 

E  -  E  E  E  1  ■* 

>■■  .  ,  1  s  in  r'.  ; 

p  q  j=--(p.q)+^  0  > 


s  i 


>: . 

2  f  2 

V.  /  C  si 

°J 


sin  (k  .)]  d>: '  + 

c  J-l 


M-l 
>:j+l 


l 


Crt.:VVl-  »•))  ««• 


M(p,q)-1 

p  q  j=N(p.q)+i 


j  pq  .  K ... 

j  sink( 


x 


j+1' 


-  2k  cos  k'  (Cd  -  Gd) |  ,  + 

x  x  1  x  -  y 

j 


f  (z  '  -z  ) 

(3 .2^ 

The  location  of  the  voltage  generator  determine?  which  of  these  tv- 
expressions  should  be  used  for  the  formation  of  the  matrix  equation. 

For  the  case  of  Figure  2,  where  the  excitation  is  or.  an  embedded  dipole, 
equation  (3.24)  is  substituted  into  (3.7).  By  using  equations  (3.3’'. 
(3.15)  and  (3.22)  one  can  determine  in  a  st ra igb.t f orvard  manner  that  tht 
elements  of  the  generalized  impedance  matrix  are  giver,  b 


+  k(Gd-Cd):  ,  +k(Gd-C,d)’  , 

X  X  «X  .  -  X  *«x . 

J-l  J+1 


'rinV-r'2  £;--h> 


E 

E 

r^k+a+l  f 

,  "  / 

a-  -1,0 

a'-  -1,0  * 

*  X  .  /  >. 

k+u  ^  ; 

d>:  * 
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4 


sin  {k|£xa  +  (xk+a+1  -x)  | } ■ 


sin  {k|txa*  +  (xj+a+1  -x')|) 


)2  +  (y-y')2  ;  z;  z’j 


f  k+a+1 

+  >  y  k  U-(2  cosk£x+l)  *£(3’)]  •  I 

a- -1,0  a'— 1,0,1  J\+a 


dx 


sin  { k [ t  a  +  (x.  ,  -x) !  } 

x  k+a+1 


(/;-x- 


)~  +  (>•->•’ )‘  ;  z;  z' 


-Gd  (  /T>:->:.+a,>2  +  (y-y*)2  ;  z;  z’) 

Thi  matrix  equation  then  takes  the  final  forr 


z’rv/pq 

ipq 

kj 

k>: 

ev 

kj 


,wv 


(3.25) 


(3.26) 


where  V,  .  are  the  elements  of  the  excitation  vector, 
kj 

For  the  excitation  of  the  driven  dipole  shown  in  Figure  3.2  one 
can  use  a  voltage  gap  generator  or  a  frill  of  magnetic  current 
[36]  located  at  the  feed  point.  From  these  two  ways  the  second  one 
although  it  has  an  exact  (or  nearly  so)  model  of  the  physical  gap  at 
the  feed  point,  it  complicates  the  problem  especially  in  the  case  of 
a  printed  dipole  since  half  of  the  dipole's  cross  section  is  in  the 
air  and  half  embedded  ir.  t  ne  dielectric.  For  this  reason  a  voltage  gap 
generator  is  considered  at  the  feed  print  and  all  the  elements  of  the 
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EVALUATION  OF  JI^E  SOMblER  FK  LD  _TY_PE  INTEGRALS 
4 - 1 .  SINGULAR  POINTS  AND  RELATED  SURFACE  WAVES 

I 

As  shown  in  Chapter  3,  the  elements  of  the  generalized  impedance 
matrix  consist  of  a  summation  of  double  and  triple  integrals.  One  of 

these  integrals  is  a  semi-infinite  integration  resulting  from  the  i 

Green's  function.  The  integrand  of  this  integral  is  an  expression  of  j 

transcendental  functions  multiplied  by  Bessel  function  of  the  first-  j 

kind  and  zero-order.  Furthermore  this  is  a  Sonmerf eld-tvpe  integral  and  j 

the  existence  of  essential  singularities  in  its  integrand  necessitates  j 

j 

very  careful  treatment.  In  this  chapter,  the  computation  of  these  inte- 

i 

grals  will  be  shown  explicitly.  In  addition,  the  approximations  employed 
will  be  justified  and  an  estimation  of  the  error  intoduced  by  these 
approximations  will  be  made. 

From  equation  (3.25)  it  is  obvious  that  all  the  singularities  of 
the  integrals  considered  here  result  from  the  Green's  function.  Further¬ 
more  as  was  shown  in  Chapter  2,  the  integrand  in  the  se~  i-inf  in.it  c 

of  the  parameter  >  through  the  radicals 

(4.1' 

(4.2'' 

functions  of  the  complex  variable  .  However, 
u  does  not  affect  the  sing! e-valuedness  of  the 
involving  the  radical  u  art.  even  functions  cf 
brand,  cut  contributed  h-  the  radical  u 

c 

of  the  branch  cut  is  not  arbitrary.  Its 


integrals  is  a  function 

*>  *»  • 
u  =  -  rr 


which  are  double-valued 
the  sign,  of  the  radical 
integrals,  as  the  terms 
.  Therer'crt,  only  the 
considered.  The  choice 


2  b 


is 


position  is  determined  by  s ingl c-val ui-dness  and  convergence  of  the 
integrals,  as  well  as  the  outgoing  wave  character  of  propagation. 

At  first,  by  considering  the  requirement  of  an  outgoing  wave  it 
can  be  shown  (see  Appendix  D)  that  the  real  and  imaginary  parts  of  the 
integration  parameter  X  must  satisfy  the  following  relations 

Re [ X ]  >  0  (4.3) 

Im[X]  >0  (4  .4) 

A  possible  position  of  the  branch  cuts  governed  by  these  two 
relations  is  shown  in  Figure  4.1.  For  a  given  complex  X,  u  can  be 
written  as  follows: 

u  =  />Vk~  •  />+k  (4.5) 


with 


j^/2 


>  X  -V.  =  eJ '  1 

c  1 


. -  , —  v  /*> 

//+r.  =  >  r _  e3^^ 

c  2 


From  (4.6),  (4.7)  and  (4.5)  one  obtains 

r -  j(t.  +  <,)/2 

u  =  rr.r,,  e  l  2 
o  l  2 


(4.6) 

(4.7) 

(  4  .  E  ) 


Fror  Figure  4.1  it  can  be  observed  that,  for  all  '•  ir.  the  first  qua¬ 
drant,  the  angles  : ^  and  C?  are  both  positive,  with 


0  <c1<  90p  +  f 


((4.0) 


and 


qriO  _  : 


(4.101 


B%  considering  (4.9)  and  (4.101  the  inequa'-it'- 

+  i  , 


•  t 

C  i  -- 


resu: : 
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This  means  that  both  the  imaginary  part  and  the  real  part  of  the  radical 
uq  are  positive.  Thus  the  choice  of  branch  cuts  as  shown  in  Figure  4.2 
satisfies  conditions  (4.3)  and  (4.4).  This  choice  also  imposes  the 
following  relations: 

Re(u  )  >  0  Im(u  )  >  0  (4.11) 

o  o 

Re(u)  >  0  Im(u)  >  0  (4.12) 

In  Chapter  2  it  was  shown  that  all  the  components  of  the  dyadic 
Green's  function  were  of  the  form 


a0) 


N(X,b,h,cr) 

DU7b,h,er) 


dX 


(4.13) 


N(X,b,h,tr)  is  a  complex  function  with  singularities  of  order  less  than 
one  (at  most),  while  the  function  Dp.b.h.t^)  is  of  the  form 


D(X,b,h,£r)  -  f1  (X ,b)  •  f2  (X,b) 


(4.14) 


or 

DO,b,h,e  )  ■=  f1  (X  ,  b)  (4.15) 

with  f^  (X ,b) ,  f2  ( X  , b )  given  by 

f.  (X,b)  *=  u  sinh(ub)  +  u  cosh(ub)  (4.16) 

1  o 

f2  (X,b)  *  tr  u^  cosh(ub)  +  u  sinh(ub)  (4.17) 

The  integrand  in  equation  (4.13)  has  poles  whenever  the  function 
DO,b,h,e  )  becomes  zero  or  when  either  one  of  the  functions  f^  (X,b), 

f2  (X,b)  become  zero.  The  zeros  of  these  two  functions  lead  to  surface- 
wave  modes.  Particularly,  the  zeros  of  f^  (X,b)  give  rise  to  TE 
surface  waves  while  the  zeros  of  f0  (>,b)  correspond  to  TV  surface 
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waves.  Assuming  that  k,  kQ  are  real  numbers  (lossless  dielectric) 
the  roots  of  equations 

u  •  -  u  coth(ub)  (4.16) 

and 


£ruo  ■  -  u  tanh(ub)  (4.19) 

lie  in  the  strip  ko  <  Re(X)  <  k  and  are  all  real.  Graphical  solutions 
for  the  roots  of  these  two  equations  are  shown  in  Figures  4.3  and  4.4 
with  the  circles  having  radii  a  given  by 


a2  =  (*^-1)  (k  b)2  (4.20) 

r  o 

From  these  figures  it  is  obvious  that  a  surface-wave  mode  exists  only 
if  a  is  greater  than  y.  As  kQb  increases,  the  number  of  proper  and 
improper  roots  increase  continually.  The  proper  roots  characterized  by 
positive  real  values  for  uq  and  positive  imaginary  values  for  u  lie  on 
the  proper  Riemman  sheet  while  the  improper  roots,  characterized  by 
negative  real  values  for  uo  and  positive  imaginary  values  for  u,  lie 
on  the  improper  Riemman  sheet.  Because  of  their  location,  only  the 
proper  roots  yield  surface  wave  modes  and  these  have  a  significant 
effect  on  the  input  impedance  and  radiation  characteristics  of  printed 
antennas. 


4-2.  NUMERICAL  INTEGRATION  OF  THE  INTEGRALS 

The  integrals  in  the  expression  for  the  elements  of  the  generalized 
impedance  matrix  can  be  put  in  the  following  two  forms: 


and 


-  f  k+dx  *<x)  |  J+1  dx’c(x’)  f  Jc(Xc1)  f(>,b,tr)  Xd>. 

Jxk  Jx  j  Jo 


(4.21) 
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Figure  4.3:  Graphical  Solution  of  the 
Equation  ueb=[u bj  cotj ub| 
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4 


Figure  4.4:  Graphical  Solution  of  the 
Equation  cfu#b  =  Jubjtanjub  j 


with 


♦<„  r 

Jx,  Jo 


J o^02)  f(X  ,b,cr)  XdX 


l(x-x')‘  +  (y-y’T] 


P2  -  [(x-x  )2  +  (y-y')2]*5 


(A. 22) 

(4.23) 

(4.24) 


The  infinite  integration  is  performed  along  the  real  axis  and  is  com¬ 
pleted  in  two  steps: 

1)  Numerical  integration  over  the  interval  [0,A]  where  A  satisfies 
the  relationship  coth  (  fTT#  b)  •  1 

2)  Analytical  evaluation  of  the  tail  contribution  which  is 
actually  the  integration  over  the  path  [A,"  ). 

These  two  integrations  will  be  analyzed  separately  because  they 
have  to  be  treated  in  very  different  ways. 

4-3.  INTEGRATION  OVER  THE  INTERVAL  [0,A] 

The  integrals  given  by  equ.  (4.21)and  (4.22)  are  now  in  the  form 


rx j+i 

dx  $(x)  /  dx'  o(x 


rv+i 

*  I  d> 

Jxk  Jx , 

rk+1  f 

dx  $(x)  / 


>/ 


J0(>:1)  f(X,b,tr)  XdX 


Jo(Xc2)  f(X,b,e  )  XdX 


Consider  that  the  function  F(c .b.c^)  is  given  by 


F(: 


*B’er>  '  / 

J  o 


J0(Xc)  f(x,b,cr)  XdX 


(4.23) 


(4.26) 


(4.27) 


In  order  to  find  the  first-order  derivative  of  the  function  FL'.b.e^) 
with  respect  to  x,  both  sides  of  equation  (4.27)  are  differentiated  and 
the  order  of  integration  and  differentiation  is  interchanged,  yielding 


d  F( 
dx 


?,b,er)  fA  d  JQ(\P) 

X  Jo  ~d*~~ 


f  (X,b,er)  XdX 


(4.28) 


Since  P  is  the  distance  between  the  source  and  the  observation  point, 
given  by 

p  <=  t(x-x’)2  +  (y-y')2]5*  (4.29) 
the  derivative  of  this  zero-order  Bessel  function  (equation  4.28)  can 
be  written  as  follows: 


d  Jo(Xp) 
dx 


J,  (Xp) 

X  — -  (x-x') 

P 


(4.30) 


Equations  (4.30)  and  (4.28)  result  in  the  following  relation: 


d  F(c,b,cr) 
dx 


(XC)  f(X,b,E  )X  dX 


(4.31) 


From  (4.31)  it  is  obvious  that  the  derivative  of  F(p,b,tr)  with  respect 
to  x  can  take  large  values  for  specific  b  and  er  which  in  turn  in¬ 
creases  considerably  the  error  of  integration  in  (4.25)  and  (4.26). 

In  order  to  avoid  this  undesirable  integration  error  it  has  been 
found  necessary  to  interchange  the  order  of  integration  so  that  the 
space  domain  integration  can  be  performed  first.  Thus,  equations  (4.25) 
and  (4.26)  take  the  form 


rA  pk+i  f 

[j  *  I  dX  X  f(X,b,cr)  /  dx  $(x)  / 

Jo  Jx  1,  -V 


k+1  /  "J+l 

dx  $(x)  /  dx'  c(x')  J_(X0l) 


(4.32) 


j  d>  X  f(X,b,cr)  r k+k  *(x)  J_ (X : ,) 

Jo  Jx\, 


"2  ”  I  “  "  I  T'"'  "o'  "2 

fo  J*  k 

Because  of  the  way  A  has  been  defined,  A  is  much  larger  than  k. 


(4.33) 


Therefore  as  X  moves  from  0  to  A,  it  passes  through  the  values  kQ  and 
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k,  causing  the  radicals  and  u  to  alternate  between  pure  real  and  pure 
imaginary  values.  This  fact,  as  well  as  the  existence  of  poles  between 
kQ  and  k,  make  it  necessary  to  divide  the  integration  interval  into  the 
following  subintervals: 

1)  0  <  X  <  kQ :  The  integration  over  this  interval  is  performed 
numerically,  using  a  modified  Romberg-quadrature  method  [39], 
[40]  for  the  integration  with  respect  to  X  and  Gaussian- 
quadrature  with  fixed  points  [41]  for  the  double  and  single 
x  and  x ' -integrat ions .  The  integrals  contribute  to  the 
radiation  resistance  and  to  the  reactive  part  of  the  input 
impedance . 

2)  kQ  <  X  <  k:  For  the  integration  over  this  interval  a  singu¬ 
larity  extraction  technique  is  used  (see  Appendix  E)  which 
transforms  the  integral  into  a  finite  series  plus  an  integral 
of  a  slowly  varying  function.  This  finite  series  gives  the 
contribution  of  the  surface  wave  modes  and  the  number  of  its 
terms  depends  on  the  thickness  of  the  dielectric  as  well  as 

the  dielectric  constant  z  . 

r 

3)  k  <  X  <  A:  Numerical  integration  is  again  invoked  here  in 
exactly  the  same  way  as  it  is  performed  in  the  first  sub¬ 
interval  . 

4-4.  TAIL  CONTRIBUTION 

In  this  case  the  integration  with  respect  to  X  is  extended  along 
the  interval  [A,  ®) .  The  use  of  the  equality 

coth  (A2  -  k^  b)  -  1  for  X  £  A  (4.34) 
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simplifies  the  integrand  to  a  form  which  involves  radicals  of  uq,  u. 

The  order  of  integration  is  not  changed  in  the  integrals  (4.21)  -  (4.22) 
but  the  Infinite  integration  is  evaluated  analytically  by  making  some 
approximations  which  are  found  to  introduce  a  negligible  error.  From 


these  results  the  integrations  obtained  by  the  method  of  moments  are 
computed  numerically.  A  technique  of  extraction  of  the  fast  varying 
terms  is  adopted  to  obtain  better  accuracy  and  to  reduce  computer  time. 

Using  the  relation  (4.34)  the  Green's  function  as  it  appears  in 
Pocklington ' s  Integral  Equations  (2.52)  and  (2.53)  is  given  by 


G  (tail)  -  /  J  (Xo)  e"uoZ  ~ -  XdX 
X  I  o  u  +u 

J  A  ° 


G(tail)  *  I  JQ(Xc)  e-uoZ  euz 


G>11)  -  j 


'  r_i _ i_i 

u  -Hi  e  u  +u  I 
L  o  r  o  J 


z  >  0 
z'  <  0 

z  >  0 


(4.35) 


(4.36) 


/*oo 

/  J(Ap) 

-u  j  z+z ' |  u-u  -u I Z+Z '  1 1 

e  1  o  (  e  '  ' 

u  u  Mi  u  J 

1  A 

u  o  J 

z'  <  0 

XdX  z,z’  <  0  (4.37) 
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G(tail)  -  I 


G  (tail) 
x 


[  1  ~e  x  (  A)  ]  [1-€2(A>]  (  Er+D  [  1  ~e  2  ( A)  111—  Ej(A)  ] 


(A) 


2[l+Cl(A)J 


rrh 


♦  1 


h ' h 


-  -  F 
P  o 


(4.40) 

(4.41) 


G  (tail)  - 


where 


11  1  1 

1  1  F 

4  1+e^A)  2  (l+er)[l+e3(A)] 

I*!  1 

[(x-x’)2  +  (y-y')2]^ 


itA 


Fo  -  A  J0(A.  )  -  -f  [J1(Ap)  H0(Ap)  -  Jq(Ap)  H^Ac)] 


(4.42) 


(4.43) 

(4.44) 


2  ** 


\  -  [(z+z’)2  +  p2  [1+€4(A)]‘] 


R,  -  l(z-z')2  +  c2  [l+c,(A)j‘J 


2,** 


(4.45) 

(4.46) 


A 


//A‘-k‘ 

Jo  (tc >tl+e4(A)I)  e-t ! z+z ' I  dt 
o 


(4.47) 


rVA  _k  -t'z-z’* 

F 2=1  Jo  (totl+E4(A)])  e  '  '  dt 


(4.48) 


The  integrals  (4.47)  and  (4.48)  are  evaluated  numerically  while 
the  functions  c^(A) ,  e2(A),  tj(A]  and  c^(A)  are  correction  functions 
because  of  the  error  introduced  from  the  neglected  higher  order  terms 
in  the  binomial  expansion.  For  a  value  of  A  approximately  equal  to 
the  error  is  no  larger  than  10  ^  of  the  value  of  the  tail  contri- 
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button  and  it  can  be  neglected. 


By  substituting  these  expressions  for  the  tail  part  of  the  Green's 
function  into  equations  (4.21)  and  (4.22)  and  by  using  a  technique  for 
extraction  of  the  fast  varying  terras  one  can  show  that  the  integrals 
are  reduced  to  a  summation  of  integrals.  The  integrads  are  summations 
of  slowly  varying  functions  and  a  finite  series  of  logarithmic  func¬ 
tions  resulting  from  the  integration  of  the  terms  and  in 

equations  (4.39)  to  (4.42). 


» 
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CHAPTER  5 


NUMERICAL  RESULTS 

5-1 .  DESIGN  PROCEDURE  FOR  MICROSTRIP  DIPOLES 

This  chapter  presents  design  procedures  for  microstrip  dipoles 
printed  on  or  embedded  In  the  dielectric  substrate.  Numerical  compu¬ 
tations  have  been  performed  for  very  thin  wire  dipoles  with  a  radius 
a  *  0.0001  Xq  .  The  dipoles  are  center-fed  by  an  in  phase  unit  voltage 
delta  gap  generator  except  for  one  particular  case  where  the  dipole  is 
fed  asymmetrically.  All  the  dimensions  presented  are  normalized  with 
respect  to  the  free  space  wavelength  Xq.  Due  to  an  assumed  time 
dependency  of  eJwt,  inductive  reactance  is  positive  in  all  plots.  The 
material  given  here  relates  the  antenna  geometry  (dipole  length,  sub¬ 
strate  thickness,  dielectric  constant,  dipole-ground  plane  distance  and 
feed  point  location)  to  antenna  characteristics  (resonant  length,  reso¬ 
nant  resistance,  current  distribution  and  bandwidth).  The  presentation 
of  the  numerical  results  is  completed  in  three  steps;  at  first  a 
dipole  printed  on  the  dielectric  interface  is  considered  and  its 
characteristics  are  discussed  in  terms  of  the  dielectric  constant  and 
substrate  thickness.  After  that,  this  dipole  is  considered’  as  being 
embedded  in  the  dielectric  and  the  change  in  its  performance  is  studied. 
At  the  position  b'  ■  the  dipole  characteristics  are  considered  in 
terms  of  the  dielectric  constant  and  the  substrate  thickness  and  are 
compared  to  the  corresponding  values  for  the  printed  dipole. 
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5-2. 


DIPOLE  PRINTED  ON  THE  DIELECTRIC  INTERFACE 


a •  RELATIVE  DIELECTRIC  CONSTANT  VARIATION 

One  of  the  most  Important  characteristics  of  a  given  dipole  is  its 
input  impedance.  Figures  5.1,  5.2  and  5.3  show  the  real  and  imaginary 
parts  of  the  input  impedance  when  the  thickness  of  the  dielectric  sub¬ 
strate  is  equal  to  0.1016Xq.  The  relative  dielectric  constant  with 
values  of  2,  10  and  35  permits  the  excitation  of  one,  two  and  three 
surface  waves  respectively.  From  these  figures  as  well  as  additional 
data,  it  has  been  concluded  that,  as  the  relative  permitivity  increases, 
the  reactance  becomes  increasingly  capacitive.  The  indicated  trend  is 
that  fewer  resonances  occur  with  larger  e^.  In  addition  the  input 
resistance  curves  show  an  increasing  number  of  oscillations  and  a 
decreasing  maximum  value.  The  latter  observation  implies  less  and  less 
energy  radiated  into  space,  i.e.  a  decrease  in  the  radiation  efficiency 
of  the  antenna,  since  energy  is  trapped  in  the  dielectric  substrate  in 
the  form  of  surface  modes.  Figure  5.4  shows  how  the  resonant  length 
Lr  varies  as  a  function  of  er>  The  effect  on  the  resonant  length  of 

each  new  excited  surface  wave  is  expressed  as  a  discontinuity  in  the 

3Lr(er) 

derivative  — — -  which  happens  exactly  at  the  transition  points. 

r 

Figure  5.4  is  also  another  interpretation  of  the  decreasing  radiation 
efficiency  of  the  antenna  since  Lr  decreases  with  increasing  e^. 

Figure  5.5  shows  the  current  distribution  on  a  dipole  of  length  L  * 

0.65  X  with  substrate  thickness  b  *  0.1016  X  and  for  e  equal  to  2, 

10  and  35.  The  increasing  er  (from  2  to  35)  results  in  a  larger 
number  of  zeros  of  the  current  distribution  and  more  electromagnetic 
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Figure  5.3:  Input  Impedance  for  a  Printed 

Dipole  with  <r=35  and  b=0.10l6\o 
Resonant  Length  Lr  =  0.1025\c 
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Figure  5.4:  Resonant  Length  vs.  Relative 
Dielectric  Constant  cr  for  the 
Printed  Dipole  with  b=0.l0l6\ 
I:  One  Surface  Wave 
II:  Two  Surface  Waves 


III:  Three  Surface  Waves 


Red)  in  mA$ 


lm(l)  in  mAs 


Figure  5.5:  Current  Distribution  on  a 

Printed  Dipole  with  b=0. 1016V 
and  Length  L  =  0.65V 


energy  trapped  in  the  dielectric  which  again  means  less  radiated  power. 


b .  SUBSTRATE  THICKNESS  VARIATION 

As  in  section  5.2a,  Figures  5.6,  5.7  and  5.8  show  the  real  and 

imaginary  parts  of  the  input  impedance  when  the  relative  dielectric 

constant  is  equal  to  2.35  while  the  substrate  thickness  b  takes  the 

values  0.1016  X  ,  0.25  X  and  0.46  X  causing  one,  two  and  three 
o  o  o 

surface  waves  respectively  to  be  excited.  As  the  substrate  thickness 
increases,  despite  the  fact  that  more  surface  waves  are  excited,  the 
possibility  of  many  resonances  is  not  affected.  However,  the  increasing 
substrate  thickness  makes  the  resonant  length  varying  as  shown  in 
Figure  5.9.  In  the  region  of  one  surface  wave  (I)  the  resonant  length 
increases,  reaches  a  maximum  for  b  approximately  equal  to  0.32  Xq 
(region  II)  and  continues  decreasing  as  b  becomes  larger.  As  in  the 
case  where  the  dielectric  constant  changes,  here  again  the  effect  of 
each  new  excited  surface  wave  on  the  resonant  length  is  expressed 
as  a  step  discontinuity  in  the  derivative  — 5 -  .  The  definition 

°lT 

of  the  bandwidth  BW  is  assumed  to  be  given  by 


BW 


2  1 


(5.1) 


where  euj,  are  the  3dB  points  of  the  response  of  the  dipole  considered 

by  its  equivalent  circuit  (see  Appendix  F) .  The  bandwidth  of  the  dipole 

and  its  resonant  resistance  as  functions  of  the  dielectric  constant  are 

shown  in  Figures  5.11  and  5.12.  In  these  two  figures,  the  fact  that 

the  bandwidth  curve  follows  in  shape  the  resonant-length  curve  is  quite 

interesting.  Both  the  resonant  resistance  and  the  bandwidth  have  their 

minimum  value  at  b  *  0.38X  and  their  maximum  value  at  b  *  0.22X  while 

o  o 
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Figure  5.6:  input  impedance  for  a  Printed 

Dipole  with  cr  =2.35  and  b=0.10l6\e 
Resonant  Length  Lr=0.3575\o 


Figure  5.8:  Input  Impedance  for 
Dipole  w 
Resonan 


b( 

Figure  5.9:  Resonant  Length  vs.  Substrate 
Thickness  b  for  a  Printed 
Dipole  with  cr  =  2 . 3 5 
I:  One  Surface  Wave 
II:  Two  Surface  Waves 
III:  Three  Surface  Waves 


Figure  5.10:  Current  Dietribution  on  a 

Printed  Dipole  with  cr=2.35 
and  Length  L*0.65\o 


(ohms) 


Figure  5.11:  Resonant  Resistance  vs.  Substrate 
Thickness  b  for  a  Printed  Dipole 
with  cr  =  2 .3 5 


the  maximum  resonant  length  occurs  between  these  two  values.  The 
current  distributions  on  a  dipole  of  length  L  ■  0.65  Xq  with  “  2.35 
and  for  b  equal  to  0.1016  Xo,  u.25  Xq  and  0.46  Xq  are  given  in  Figure 
5.10.  From  these  figures  it  can  be  observed  that  the  increasing  number 
of  surface  waves  coming  from  changes  in  the  value  of  b  does  not  affect 
the  current  distribution  so  much  as  it  does  when  er  changes.  Here  again 
larger  substrate  thickness  means  more  electromagnetic  energy  trapped  in 
the  dielectric  and  less  radiated. 

5-3  .  DIPOLE  EMBEDDED  IN  THE  SUBSTRATE 

a .  VARIATION  OF  DI P0LE_  EMB EDDING  POSITION 

Figure  5.13  shows  the  real  and  imaginary  parts  of  the  input  impe¬ 
dance  of  a  dipole  printed  on  the  dielectric  interface  with  relative 

dielectric  constant  e  ■  3.25  and  substrate  thickness  b  ■  0.1016  X  . 

r  o 

If  the  dipole  is  embedded  in  the  dielectric  substrate,  the  input  impe¬ 
dance  will  change  as  shown  in  Figures  5.14,  5.15  and  5.16.  From  these 
figures  it  can  be  observed  easily  that  the  shape  of  the  curves  remains 
the  same  while  the  maxima  and  minima  are  shifted  to  smaller  values  of 
the  dipole  length.  Figure  5.17  shows  the  resonant  length  as  a 
function  of  the  dipole-ground  plane  distance.  The  change  in  the 
current  distribution  as  a  dipole  of  length  L  *  0.65  Xq  is  placed  closer 
to  the  perfect  conductor  plane  is  shown  in  Figures  5.20  and  5.21  for 
four  different  cases:  b’  *  b,  b'  ■  b,  b'  ■  ^  b  and  b1  «  |  b 

where  b*  is  the  dipole-ground  plane  distance.  Also  Figures  5.18  and 
5.19  show  the  resonant  resistance  and  the  bandwidth  BW  as  functions 
of  b-b',  the  distance  of  the  dipole  from  the  air -dielectric  interface. 


Figure  5.18:  Resonant  Resistance  vs.  Embedding 
Distance  b-b'  for  a  Dipole  with 
cr=3 . 2 5  and  b=0.10  16\o 


6  ] 


Figure  5.19:  Bandwidth  vs.  Embedding  Distance 
b-b'  for  a  Dipole  with  cr=3.2  5 
and  b=0.1016\o 


Red)  In  mAi  lm(l)  In  mAi 


Figure  6.20:  Current  Distribution  for  a  Dipole 
Embedded  in  the  Substrate  with 
cr=3. 2 5 ,b-0 . 1  016  and  Length  L  =0.65\ 

(a) :  Embedding  Distance  b-b'=0 

(b) :  Embedding  Distance  b-b'=b/l00 
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Red)  In  mAs 


ImO)  In  mAi 


> 

_ I _ _ 

4  2  0  2  4 


Figure  5.21:  Current  Distribution  for  a  Dipole 
Embedded  In  the  Substrate  with 
cr  =3.25, b=0. 1016  and  Length  L=0.65\o 

(c) :  Embedding  Distance  b-b'=b/10 

(d) :  Embedding  Distance  b-b'=b/3 
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Iiu  high  slope  of  those'  two  curves  as  well  as  of  the  resonant  length 
curve  (Figure  5.17)  for  small  values  of  b-b'  is  characterf ~t ic  of  the 
air-d ielectr ic  substrate  electrical  discontinuity, 
b .  REIATIVE  DIELECTRIC  CONSTANT  VARIATION 

We  consider  now  the  dipole  shown  in  Figure  5.22.  The  thickness  of 

the  dielectric  substrate  is  equal  to  0.1016  >  and  when  it  is  combined 

o 

with  e^  =  2  and  «  10  it  permits  the  excitation  of  one  and  two 
surface  waves  respectively.  The  real  and  imaginary  parts  of  the 
input  impedance  for  these  two  cases  are  shown  in  Figures  5.23  and  5.24. 
From  these  figures  it  can  be  observed  that  for  the  case  of  two  surface 
waves,  as  well  as  for  the  case  of  one  surface  wave,  many  resonances  are 
possible  as  the  length  of  the  dipole  increases. 


This  does  not  happen  when  the  same  dipole  is  printed  on  the  di¬ 
electric  interface.  Figure  5.25  shows  how  the  resonant  length 
varies  as  a  function  of  c  for  three  different  positions  of  the  dipole; 

a)  printed  on  the  interfaced),  b)  embedded  in  the  substrate  at  a 
2  B 

distance  -j  E  (2), and  c)  y  (3) from  the  ground  plane.  It  is  quite 
interesting  to  note  that  as  the  dipole  enters  the  dielectric  substrate, 
its  resonant  Length-Dielectric  Constant  curve  does  not  change  shape 
but  moves  to  smaller  values  for  the  resonant  length.  This  means  that 
curve  5.17  maintains  its  shape  as  the  dielectric  constant  1 

r 

changes.  For  this  reason  it  is  expected  again  that  the  effect  of  each 


new  excited  surface  wave  on  the  resonant  length  will  be  expressed  as 

?Mtr) 

a  discontinuity  in  the  derivative  — r -  and  the  curve  L  *  f(e  )  for 

c  Z  TV 

r 

fr  larger  than  10  will  exhibit  a  similar  variation  as  the  one  for  the 


printed  dipole  (Figure  5.5).  Figure  5 . 2 f  shows  the  current  distribution 
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Figure  6.22:  Dipole  Embedded  in  the  Dielectric 
Substrate 


U 

L _ _ _ 


.6 


L(X.) 

t  Impedance  for  a  Dipole 
edded  in  the  Substrate  with 
'2,cr  =  2  and  b=0. 1016V, 


Figure  6.25:  Resonant  Length  vs.  Relative  Dielectric 
Constant  cf  for  a  Dipole  with  b=0.1016\ 
o — o  Dipole  on  the  Interface 
e— e  Dipole  in  the  Substrate  b'=b  1 2 
4 — *■  Dipole  In  the  Substrate  b'=2 b / 3 


-005 


on  tin.  dipole  of  Figure  5.22  with  its  length  equal  to  0.65  and  for 

two  values  of  e  (e  *  2,  10).  The  increase  in  the  number  of  zeros 
r  r 

resulting  from  the  increasing  relative  dielectric  constant  is,  as  in 
the  case  of  the  printed  dipole,  characteristic  of  the  substrate  proper- 
t  ies . 


c.  Sl’BSTRATE  THICKKESS  VARIATION 

The  dipole  of  Figure  5.22  will  now  be  considered  again.  For  this 

dipole  with  a  relative  dielectric  constant  for  the  substrate  equal  to 

2.35,  the  real  and  imaginary  parts  of  the  input  impedance  for  three 

different  values  of  the  substrate  thickness  are  given  in  Figures  5.27, 

5.28  and  5.29.  These  three  values  of  b  (0.1016  >  ,  0.25  )  and  0.46  >  ) 

o  o  o 

cause  one,  two  and  three  surface  waves  to  be  excited.  As  shown  in 
Figure  5.30,  the  resonant  length  as  a  function  of  the  substrate  thick¬ 
ness  goes  through  a  minimum  in  region  1  (one  surface  wave)  and  continue? 

increasing  in  region?  II  (two  surface  waves)  and  III  (three  surface 

8Lr(b) 

waves)  with  a  step  d iscont inuitv  in  its  derivative  -r.  —  at  the 

c  0 

point  where  a  third  surface  wave  is  excited.  Assuring  the  same  defi¬ 
nition  of  bandwidth,  as  in  section  5.2,  one  can  deter-ir.c  its  variation 
with  respect  to  substrate  thickness  b.  Ir,  addition  the  resonant  re¬ 
sistance  as  a  function  of  the  same  parameter  is  shown,  ir.  Figure;  5.31 
and  5.32.  From,  these  figures  one  can  see  that  the  resonant  resistance 


is  monotonical ly  increasing  as  h  increases  with  step  d i scont inu it ies 
?R_  fh) 

in  the  derivative  — - —  at  the  points  where  a  new  surface  wave  is 

'a  V 

excited.  The  bandwidth  B"  follows  in  shape  thi  variation  of  Kr  but  it 
exhibit?  murh  larger  slope.  Figure  5.??  show?  the  current  d i st r ibu t io^ 


o f  the  d 


ipcle  of  Figure  5.22  when.  1  *  0 . 6  ^  ^ 


=  2.35  and  b 
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Fiflure  5.30:  Resonant  Length  vs.  Substrate  Thickness 
b  for  a  Dipole  Embedde  in  the  Substrate 
with  «r  =  2.3S 
I:  One  Surface  Wave 
II:  Two  Surface  Waves 
III:  Three  Surface  Waves 


Figure  5.32:  Bandwidth  BW  vs.  Substrate  Thickness 

b  for  a  Dipole  Embedded  in  the  Substrate 
with  b'=b/2  and  tr=2. 35 


0.1016  )  ,  0.25  X  and  0.46  X  . 

o  c  o 


Figures  5.34  and  5.35  show  the  current  distribution  as  a  function 
of  the  length  of  the  dipole  (Figure  5.22)  when  b'  *  2b/3,  b  *=  0.1016  X 
and  c  -  2.35.  The  feed  point  is  not  located  at  the  middle  of  the 
dipole  contrary  to  all  the  cases  which  were  studied  until  now.  Thc- 
distance  of  the  feed  point  from  the  nearest  end  of  the  dipole  is  1/4 
of  the  dipole  length. 
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INTEGRAL  RF_P  RESEN  TAT  I ON  OF  ITiF  PRIMARY 
SOLUTION 

As  mentioned  in  Chapter  2,  the  primary  solution  is  the  particular 
solution  to  the  wave  equation  given  by: 


p>: 


jwp  -Jkr 
J  o  e  J 


Ark* 


(  A  .  1 ) 


It  is  desired  to  bring  the  representation  (A.l)  of  into  the  forr 

of  superposition  of  e  igenf unc t ions  Since  the  cylindrical  polar 

coordinates  r,  i,  z  are  involved,  these  eigenf  unc  t  ions  will  be  ir.  tf.c 
form  [31 

(a.:' 


T  /,  ,  in:  ii/k'-'“ 
,  =  J  ( >■ : )  e  e 

r.  n 


•>  h 

where  ;  =  [(y.-x’)-  +  (y-y  )]*  and  )  are  eigenvalues  which  for  the  case 
of  an  infinite  cylinder  with  finite  radius  are  restricted  to  a  dis¬ 
crete  spectrum.  However,  if  f  •*  «  (as  is  the  case  here',  '•  has  a 
continuous  spectrum.  0  <  )  <  «  corresponding  tc  the  unlimited  reriur  . 
Also,  from  the  eigenfunctions  vi  only  those  are  used  whirl,  are  inde¬ 
pendent  of  i  ;  therefore  n  must  be  zero.  From,  ti  ic  discussion  it  is 
concluded  that  the  particular  solution  (A.H  may  he  put  in  the  forr: 


d  J'  •  c 

:  I: 


J  (>..1  e~u  74V  d- 

o 


vat  r. 


u  ■=  [>‘-k“]  • 


(A. 3) 


(A. A' 


;:‘4(Z+n* :  ■ 


8  3 


r 
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is  tin  Four icr-Bessel  Transform  of  the 


In  equation  (A. 3)  A()c 
o'jkr 


-u ; ?+L 


function  -- 


and  therefore  is  given  by 


A(>.)e"ulZ+h! 


}■* 


JD(Xf)  dp 


(A. 5) 


Equation  (A. 5)  is  easily  simplified  by  considering  z  *  -  h  which  gives 


A(X) 


/. 


e‘jkf  J  (XP)  dp 
o 


In  (A. 6)  J  (Xc)  is  replaced  by  its  integral  representation 


J  (X 
o 


-•  / 


JXpcosw  dw 


and  the  relation 


AO)  = 


±  j  dv  d:  . 

i  r  dv _ 

■  -  >  3  J  -t«COSV 

results.  The  regaining  integration  in  (A.E)  gives 


(A. 6) 


(A. 7) 


(A.E) 


AO) 

If  equation  (A. 9)  is  substituted  into  (A. 3)  the  following  relation  for 


G  is  obtained 

P>: 


3  .  _  hi*  f 

4-k2  J 


J  (>.-)  e-u:z+h  id'- 

0  u 


(A. 10) 


8t 


* 


Append  1>:  E 


BERT  V  ATI  OK  OF  TRF.  COMPOKF.KTS  G  and  G 


During  the  formulation  of  Pocklington's  integral  equation  the 
following  two  relations 


3G 
_ z 

3z 


3G 

c)X 


(B.l) 


3G 


cZ 


MI 

3>: 


were  obtained.  Bv  observing  that  one  can  write  these 

cx  cx 


equations  as  follows: 


cC 


dG  =  -rfdv- 

cZ 


,  5C-d 

dC-d  *  dx’ 

oz 


or 


G  = 


/oG 

r  3cd 


Equation  (2.3-)  in  Chapter  2  yields,  for  h  =  (’ 


2-1 


(1-c 


f)  cos tj* 


m  /  .  \  “Vi  7 

Jj(^)  e  o 


(B.2) 


two 


(B.3) 


(B.4) 


(1.5) 


(E.O 


girth (uh_)  cosvTuv  ) ,  2  . . 

V  ° 


Ey  differentiating  (B.7)  with  respect  tr  r  and  by  changing  the  order 
of  integration  and  different  iatio:  one  find.-  the  following  relation: 
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3C  ju; 


~2  (l-cr)  cos 


■  r 


-uo)  Jj(Xr)  e~uoz 


sinh  (ub) _ cosh(ub) _ .  2  . 

"  fjCx.b)  •  f2(x,b) 


r* 

- 7;  (l-c)  cos<  /  u 

27ik2  r  Jo 


d  j  (Xr) 

o  -u  z 

0  d  ( >  r )  e  0 


sinh(uM  cosh(ub), 2 


Fron  the  relations 


i  =  Kx-x’)2  +  (y-y*)2]’* 


it  is  deterrined  tha: 


d3-I>  I  d  Jc(>- }  dx; 
d(V>  =  d>: '  d 


1  d  1 


A  substitution  of  (B.12)  into  (3.E)  vields  the  follovinc  relation 


V.  f"  d  J  (•  ' 

— (t  -1>  /  —4 -  e“V" 

2-i-  r  J 

C  ^  £' 

sinVu  1  ■  c  0  <■  !,iu! 


Appendix  C_ 


UNIQUENESS  AND  EXISTENCE  OF  THF  SOLUTION  OF  THE 

FUNCTIONAL  EQUATION  L  (J)*(E  x) 

- - -  op  x 


A.  Uniqueness 

Let  it  be  assumed  that,  for  the  equation  L0p( J)“(Exx) ,  two 


vectors  and  exist  such  that  [42] 
L  (J  )  -  (Ex)  and  L  (J  )  -  (Evx) 

op  1  X  op  l  X 


(C.l) 


A  subtraction  of  one  from  the  other  yields  the  conclusion  that 

L  (J.-J-)  »  0  or  L  (d)  «=  ($  (C.2) 

op  i  z  op 


where  d  is  the  solution  to  the  homogenous  equation  Lo^(j)  =  6. 
Because  of  the  physics  of  the  problem,  with  no  excitation  the 
induced  current  on  the  dipole  is  identically  zero  and  there- 

— V  -*  — ►  — V 

fore  d  »  0  which  gives  J,  =  J0.  Therefore  the  solution  is 

1  i. 

unique . 

B.  Existence 

The  homogeneous  adjoint  equation 


La  (Ja)  =  0 
op 


is  considered  with  J  in  the  adjoint  domain 
definition  of  the  adjoint  operator 


<  L  (Jl,  J*  > 


J,  L  (J“) 
op 


(  C .  3) 
By  the 

(C.4) 


or 


'  <V>. 


J,  c 


(C.51 
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(C.6) 


Equation  (C.5)  reduces  to 

<  (E  x)  ,  Ja  >  ■  0 
x 


Since  it  has  been  proved  that  the  operator  L  is  6elf 
adjoint , 


op 


op 


(C.7) 


results . 

Because  of  (C.7),  equation  (C.3)  reduces  to 

L  (Ja)  «  0  (C .8) 

op 

vhich  vas  shown  earlier  to  have  a  solution  identically  equal 
to  zero.  Therefore  equation  (C.6)  is  satisfied,  and  the 
existence  of  the  solution  of  the  original  functional  equation 
is  concluded. 


Appendix  _D 


CONDITION’S  IMPOSED  BY  THE  OUTGOING-PROPAGATION- 
CHARACTER  REQUIREMENT 

As  shown  in  Chapter  2,  both  of  the  components  of  the  Green's 
Function  in  air  include  a  factor  e  Uoz  in  their  integrands.  This 
factor,  together  with  the  exponential  time  dependance  factor  takes 
the  form 

c(z,t)  -  e"uoz  +  j‘t  (D.l) 

The  transformation 

u  *  j  k  cosv  (vc O  (D.2) 

o  J  o 

leads  to  the  following  relationship  for  > 

>  *  k  sinv  (D.3) 

c 

Assuming  that  w  *  t  +  jv  one  can  show  that  equations  CD. 2)  and  (D.3) 
become 


*  V  sir."  sinh 
o 

+  j 

k  cos" 
0 

coshv 

(D.4) 

*  k  sin:  coshv 

0 

+  j 

k  cost 
0 

sinhv 

(D.5) 

Because  of  (B.4),  equat  ion  (D.l)  car.  be  put  in  the  form. 

,  ,  Ifut  -  k  z  cost  cos-  coshv •  -k  z  sin*  sinhv  (D.f) 

c(z  ,t)  *  eJ  o  to 

In  order  tc  satisfy  the  outgoing  character  of  the  propagating  waves, 

the  following  relations  must  be  satisfied: 

sir.  t  cir>  ■  (  (D.7) 

cos*  cos':..'  (  ( P .  t ) 
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# 


From  (D.7),  (D.8),  and  from  the  fact  that  the  contour  of  integration 
with  respect  to  X  is  extended  along  the  positive  real  axis  it  is  obvious 
that 


i 


sint  >  0 

(D.9) 

sinhv  >  0 

(D.10) 

cos  t  >  0 

(D  .11) 

A  consideration  of  inequalities  (D.9) 

-  (D.ll)  into  (D.5)  indicates 

that  the  real  and  imaginary  parts  of 

the  parameter  X  satisfy  the 

relations 

Re  [>.]  >0 

(D.12) 

Im  [>.]  >  0 

(D.13) 
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Append  i>:  F. 


TECHNIQUE  FOR  THE  EXTRACTION  OF  THE  SINGULARITIES 


As  mentioned  in  Chapter  4,  the  integration  along  the  interval 

(k  ,k]  faces  the  existence  of  a  finite  set  of  discrete  essential  sin- 
o 

gularities  which  with  a  surface  wave  character,  affect  considerably 
the  input  impedance  and  radiation  characteristics. 

If  it  is  assumed  that  S  is  the  set  of  these  singularities  then 
S  is  given  by 


S  =  {x^/(x..=  root  of  f^().  ,b)  *  0]  U  {x^=  root  of  f,(>.  ,b)  *  0] 
and  i  =  1 ,  2 ,  ...  N; 


(E.l) 


A  partition  P  [t  j]  of  the  interval  [kQ,k]  is  considered  such  that 


{t  *  k  ,  t,  <  x,  ,  . . .  t  <  x  ,  t  *  k! 
o  o  1  1  n  n  n+1 


(E.2) 


The  integrals  (4.32)  and  (4.33),  with  the  > -integrat ion  extended  fror 


k  to  k,  can  be  put  in  the  fort 
o 


-t 

■A  l 


d  '  fv  *  •  » '  r  ^  rV  k  ^ 


EJ  (>:,.) 
o  k.i 


FJ  0  .  .)  = 
c  ki 


r"”-  r 

k  1 

Pk>l  , 

I  dx  •  (x^  Jcf'-2 


j+i 

dx  c(x  )  J  (•..■« 
o  1 


(E.3) 


(E  .4'\ 


(E .  5) 


Ecuatior.  fJ.3)  car  alsr  he  writter  as  fellow- 


r.  r-t-]  iV:(  -  .h.'  1 

•L  I  d  -r  — r-  - 

r  • 0  A 


(E.ei 


L 


with  <  (X,b,e  )  «=  X(X-x  )  f  (X  ,b,c  )  FJ  <Xp  ) 
r  r  r  r  o  kj 


(E.7) 


♦k^(x  ,b,c  ) 

In  (E.6)  the  quantity  —  y  is  added  and  substracted  giving 


£  / 


K  rl r+1  4kj(X,bte  )-  «*3(x  ,b,e  ) 

r  /  ax— - - — - — 1 — —  + 


kj, 


X  -  x_ 


r«0  It  X  -  x 

J  r  r 


or 


V'  f  r+1  ‘  r  2  (>  *b*  Cr)  *  •rJ(xr,b,cr) 

I  -  Z-  /  — - 1 - 5 1 - —  + 

r=0  J t  X  -  xf 

Jf*  tr+I 
• 

r 

From  equation  (E.9)  the  integral  can  be  put  ir.  the  final  forr. 


(E.8) 


(E.9) 


T  _  £  f  ^1<r-1<?’b'cr)  -  *rJ(Vb’cr^ 

r»0  /  X  -  x 

Jtr 

t  ?<*,*.<,•>  [’-  !¥-K;:-  J-  ] 

r=f  *■  T  T  J 


In  the  equation  above,  the  integrands 


(E.10) 
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are  slowly  varying  functions  In  the  intervals  (tr>  tf+1l  and  the  inte 
grals  can  be  evaluated  using  a  Gaussian-quadrature  integration  method 
with  four  fixed  points  {36], 


Append  ix_  _F 

BAN DW IDTH  OF  A  _DIPOLE  PRINTEPJ3N  OR 
EMBEJDDEP_IN  THE JMELECIRiC. _SJIP-SIMI£ 


We  consider  the  dipole  by  its  equivalent  circuit  (Fig.  F-l) . 


The  impedance  seen  by  the  voltage  generator  is  given  by: 


Z.  =  R  +  j  (wL  -  —  ) 
in  J  wC 


(F.l) 


With  X  being  the  imaginary  part  of  the  input  impedance,  equation  (F.l) 
can  be  written  as  follows: 


with 


Z.n  =  R  +  j  X 


>:  «=  lL  - 

L  C 


(F .  2) 


A  differentiation  of  both  sides  of  equation  (F.3)  with  respect  to 
frequency  yields  the  relation 


T  ‘  L 


t1  +  ^’] 


with  the  resonant  frecuer.cv  of  the  dipole ,  give:  r; 

r 


(F.41 


1  r  s'! 


A/. 


With  u  =  u  in  equation  (F.4),  the  derivative  of  the  reactance  with 
respect  to  frecuencv  at  resonance  is  given  by 


-c.u 

dx’ 

d.  =  •• 


(F.O 


If  .]  and  . 0  (^2  >  uT)  are  the  3  dF  point?  or.  the  response  of  the 
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toned  circuit  (Fig.  F-l)  the  quality  factor  Q  satisfies  both  the 


ret  at  ions 


and 


w  L 


(F.7) 


(F  .8) 


From  equations  (F.7)  and  (F.8)  the  relation 


U2  "  L1  =  L 


(F.9) 


results . 

If  one  substitutes  (F.6)  into  (F.9)  and  if  the  following  bandwidth 
(BV)  definition  is  assumed 


BV 


**2  -L1 


(F.1C) 


then  the  relation 


r 


is  concluded  wit'. 


(F.ll) 


d>:‘ 

d.  V 


L. 


( f . i 


where  c  is  the  velocity  of  light  in  the  medium  surrounding  the  dipole 
and  I  is  the  lencth  of  the  dipole.  Fror  equations  (F.ll)  and  (F.12'' 
it  can  be  founJ  that  the  ba^dvidtl  is  giver,  by 


BV  = 


1 

r 


(F.13) 


where  1  is  the  Resonant  Length  normal ized  with  respect  to  free  space 
wavelength. 


